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^ \ We present a crystal field theory of transition metal impurities in semiconductors 

c/j ■ in a trigonally distorted tetrahedral coordination. We develop a perturbative scheme 

d . to treat covalency effects within the weak ligand field case (Coulomb interaction dom- 
inates over one-particle splitting) and apply it to ZnO:Co^"^ (Sd^). Using the large 

Ch ' value of the charge transfer energy Apd compared to the p-d hoppings, we perform 

^ ■ a canonical transformation which eliminates the coupling with ligands to first or- 

CS| . der. As a result, we obtain an effective single-ion Hamiltonian, where the influence 
> 

CN ' of the ligands is reduced to the one-particle 'crystal field' acting on d- like- functions. 
(N ■ 

£_2 ' This derivation allows to elucidate the microscopic origin of various 'crystal field' pa- 



rameters and covalency reduction factors which are usually used empirically for the 
interpretation of optical and electron spin resonance experiments. The connection 
of these parameters with the geometry of the local environment becomes transpar- 



i-rt ■ ent. The experimentally known (/-values and the zero-field splitting 2D are very well 

O . reproduced by the exact diagonalization of the effective single-ion Hamiltonian with 



only one adjustable parameter Apd- Alternatively to the numerical diagonalization 



^ ■ we use perturbation theory in the weak field scheme (Coulomb interaction ^ cubic 

splitting S> trigonal splitting and spin-orbit coupling) to derive compact analytical 
expressions for the spin-Hamiltonian parameters that reproduce the result of exact 
diagonalization within 20% of accuracy. 

PACS numbers: 75.50.Pp,75.10.Dg,76.30.Fc 



I. INTRODUCTION 

The A^^B^^ and A^^^B^ compounds with transition metal impurities are called diluted 
magnetic semiconductors (DMS). These systems are interesting both from practical and 
fundamental points of view. Recently, they have attracted great interest as potential mate- 
rials for spintronics, the current trend of electronics which manipulates also the spin of the 
carrier, not only its charge. The challenge to physicists is the highly correlated subsystem 
of transition metal ions (TMI) demonstrating peculiar interplay of the covalency and of the 
localized nature of rf-electrons. 

Usually, the results of optical and electron spin resonance (ESR) experiments are phe- 
nomenologically describecli>2ii in terms of crystal field (CF)^'^'^ theory. It considers the 
one-site Hamiltonian for the (i-shell of TMI, where Coulomb interaction and spin-orbit (SO) 
coupling terms are added to one particle CF term of unspecified nature. The CF interpreta- 
tion reflects the symmetry of TMI environment, but the values of energetic parameters are 
taken from experiment, and, thus, say nothing about physical interactions behind them. It 
concerns not only the one-particle CF terms: phenomenological reduction factors for Racah's 
Coulomb parameters B and C, for the orbital angular momentum, and for the SO coupling 
are always introduced. 

Within the phenomenological CF theory, the hierarchy of interactions is well established 
in most complexes containing the transition metal ion surrounded by ligands. For commonly 
met octahedral and tetrahedral environments, the CF terms are divided into a large cubic 
part and much smaller low symmetry terms. For the 3d ions the SO term is one or two 
orders of magnitude smaller than the cubic terms. To decide the order of the interaction 
strength between Coulomb forces and cubic CF is much more difficult. One distinguishes 
the 'strong' and 'weak' CF cases. ^ In the 'strong' field approach,— the cubic Hamiltonian is 
diagonalized first, and its eigenfunctions serve as a basis for the subsequent consideration of 
Coulomb and remaining terms. In the 'weak' field approach, the many-body eigenfunctions 
of the single-ion Hamiltonian are exposed to the action of CF. 

The Hartree-Fock molecular orbital theory enables to take into account the covalency 
within the 'strong' field approach.^ Then, in principle, the calculation of CF parameters be- 
comes possible, the origin of the reduction factors becomes clear and they can be estimated. 
We mention here the simple analytic approach proposed in Ref. H for the CF parameters 



calculation. It is based on Harrison's parametrization of the electronic structure of solids^ 
and it enables to relate the CF parameters with the structure of the TMI environment. 

In the 'weak' field case, the account of covalency involves the Heitler- London configuration 
interaction (CI) approach.^ It is much more complicated as it works with the enlarged Hilbert 
space of many-body functions. On this way, the description of the DMS on the energy scale 
of the cubic splitting (~ 0.5 eV) was achieved,— but more fine properties (such as magnetic 
anisotropy) were not obtained, and the relation of CI approach with the phenomenological 
CF theory was not established. 

The A^^B^^ semiconductors crystallize in zincblende or wurtzite structures. In both 
structures, the TMI substituting for the A site has similar tetrahedral coordination. The A 
site has cubic point symmetry in zincblende structure, and trigonal symmetry in wurtzite 
structure. In CF theory, the cubic field of the zincblende structure is described by one 
parameter A, and the trigonal field (wurtzite structure) by three parameters A, v, v'. The 
deviation from cubic symmetry in wurtzite structure, though small (A ^ v,v'), has crucial 
influence on the magnetic properties of TMI. The Co^+ ion in trigonally distorted environ- 
ment displays a large single-ion magnetic anisotropy, which strongly affects the magnetic 
properties of Co-based DMS.^^ The anisotropy is linear in the trigonal field parameters v 
and v'. The empirical determination of v and v' is difficult and ambiguous. In optical 
transitions they are masked by a large A, and various sets give similar description,i2*i^ their 
relation with the TMI environment structure remains unclear. It is thus highly desirable to 
have a microscopic model that would at least diminish the number of adjustable parameters. 

One microscopic model of CF is well known and was the first one that appeared with 
the CF concept itself. It is the point charge model (PCM).— In this model, it is possible 
to calculate CF both in 'strong' and 'weak' field approaches, but the CF is then strongly 
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underestimated because the hybridization with ligands is neglected. The authors of Ref. 
have tried to improve the PCM and proposed to modify the TMI d-function in such a way 
that it gives the correct value of A within the PCM. In fact, they changed one set of param- 
eters (CF matrix elements) to another one (those describing their Slater- type d-function) . 
In order to explain the reduction factors, they additionally introduced phenomenologically 
the hybridization with ligands that changes the Coulomb and SO interactions but does not 
contribute to the CF splitting. 

Nevertheless, at present time, any modification of the PCM cannot be accepted as a 



physical model of the TMI in DMS. It is not adequate even for an isolated impurity, e.g. it 
cannot explain the increase of cubic splitting with the increase of covalency. Nor, it cannot 
be the starting point for studies of the interaction of TMI with the host valence band {p — d 
exchange) and between the impurities (super exchange). Both phenomena are due to the 
virtual hoppings of electrons between ligand and TMI described by the p — d Hamiltonian.i^ii 

In this work we establish the relation between the configuration interaction (CI) approach 
and the phenomenology of the CF. Below we will show that the many-body perturbation 
theory gives the possibility to consider the d-level splitting by covalency in the 'weak' field 
case at the same level of simplicity and physical transparency as for the pure electrostatic 
case (PCM). Our final goal is the derivation of the effective low energy spin-Hamiltonian, 
i.e. an effective Hamiltonian that describes the ground state manifold response to the 
applied magnetic field. We will demonstrate with an example for the Sd"^ ion in tetrahedral 
coordination that spin Hamiltonian may be derived from realistic many-body Hamiltonian 
by means of a set of canonical transformations. 

In a first step, one eliminates the coupling with ligands: i.e. we pass from realistic p — d 
Hamiltonian to CF-like Hamiltonian for d-ion only (Section II. A). We find a renormalization 
(i.e. reduction) of Coulomb and SO parameters (Section II. B). We generalize the ideas of 
Ref. |8J and give the expressions for CF parameters via d — p hopping and charge transfer 
energy difference. It is important to note that these expressions give the connection of CF 
parameters with the structure of the TMI environment (Section HI). 

Having derived the parameters for the effective CF-like Hamiltonian in Section HI, we 
calculate then in a second step the zero-field splitting 2D and the ^^-factors, i.e. the param- 
eters of the low-energy spin-Hamiltonian. The basis for the crystal field Hamiltonian of a 
3(i^ configuration including SO coupling has 120 dimensions and will be diagonalized exactly 
(details are given in the Appendix). It follows that we do not have to worry about the 'weak' 
field hypothesis for this effective Hamiltonian: indeed this procedure is equally valid if cubic 
splitting A and Coulomb interaction parameters B and C would be of the same order. 

On the other hand, in the given case of ZnO:Co, we can alternatively obtain an analytical 
closed expression which connects the parameters of the microscopic Hamiltonian with the 
parameters of the effective spin-Hamiltonian. For that purpose we use perturbation theory in 
the spirit of the 'weak' field scheme. We construct the many-body basis for a d"^ ion (Section 
IV. A) and make two subsequent transformations (Section IV. B). The first one takes into 



account the fact that the cubic sphtting is smaller than the remaining Coulomb interaction, 
and the second treats the trigonal coupling as a perturbation. In the final step we take into 
account spin-orbit coupling in order to obtain the spin-Hamiltonian observed in ESR. 

In the Section V, we apply our approach to the case of ZnO:Co and we compare the 
numerical diagonalization results with the perturbative formula finding reasonable agree- 
ment. We compare also our parameter set which was calculated microscopically with the 
phenomenological parameter sets derived before from optics and ESR. 

II. MICROSCOPIC FOUNDATION OF CRYSTAL FIELD THEORY 

A. Crystal field Hamiltonian 

In this section we show that starting from realistic p — d Hamiltonian it is possible to 
derive an effective single-ion Hamiltonian. The latter has the form of a classical crystal 
field Hamiltonian. The crystal field parameters acquire a clear microscopic meaning and 
the essential property of being calculable and connected with the structure of the TMI 
surroundings. The main point is the large energy separation of configurations rf" and dT'^^L 
(respectively n electrons in (i-shell of the TMI and n + 1 d-electrons plus one hole in the 
valence band), compared with the TMI-ligand hopping.^^ 

The appropriate Hamiltonian should be written in the basis of spherically symmetric 
functions. The basis should not necessarily coincide with the one for free ions. A more con- 
tracted basis is suitable for solids, e.g. that one of the FPLO method.-™ Without specifying 
it, we will assume the existence of an orthonormal basis of one-particle spherically symmetric 
functions localized on lattice sites. We assume also that it is a 'minimal' basis set, i.e. one 
radial function suffices for the description of one electronic shell. The existence of such a 
basis is not explicitly proved, but there are indirect evidences in favor of it. First, the FPLO 
method enables to explicitly construct non-orthogonal minimal basis of localized functions. 
Second, W.A. Harrison succeeded in describing the electronic structure of a huge number of 
compounds, assuming the existence of such basis, and empirically fitting the Hamiltonian 
matrix elements in this basis.-^ The Hamiltonian may be written as 

H = Hd + Hp + Tpd , (1) 

where Hd, Hp are local Hamiltonians for the TMI and ligands, respectively, Tp^ describes 
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electron hoppings between the TMI and ligands. In the superposition model,— the contri- 
butions of separate ligands are superimposed above each other. Further on, we will give a 
foundation for that rule, following the lines given in Ref. M- So, it suffices to consider the 
ligand at the point (0, 0, R). Then, the generalization to another geometry is straightforward. 
In the zero-order Hamiltonian we include the diagonal one-particle terms and dominant 
Coulomb interaction 

^0 = edNd + epNp + Ud + Up, (2) 

where 



N^^EE 



n. 



m,s ) ^m 



^m,s^'m,s ; tVi 



^{N! 



N, 



s m=-l 

ed, ep are the one-particle energies of d and p states. Ad and Ap are the corresponding Racah's 
parameters; the operator cj^^ = dj^siPms) creates an electron with the one-particle basis 
d{p) wave function with angular momentum and spin projections m and s on the TMI and 
on the ligand site respectively. In the ground state of Hq, the d-shell of the TMI contains 
n electrons and the ligand has the closed p-shell with rip = 6 electrons. The hopping 

Hamiltonian 

1 

Tpd = Z^ Z_^ tpdm \dm,sPm,s + Pm,s^rn,sJ (3) 

s m=-l 

couples configurations with different numbers of d-electrons. The most important is the 
coupling between the configurations d'"'p^ and d'^'^^p^. The hybridization with the conduction 
band depends on second nearest neighbor d—s hopping matrix element and may be neglected. 
In the two center approximation, the hopping Tpd (Eq.®) is diagonal over the angular 
momentum projection indices m and m' due to the symmetry with respect to the TMI- 
ligand axis. We perform a canonical transformation which eliminates the hopping to first 
order 



Heff = exp{-W)H exp{W) ^H + 
We choose for this the operator W in the form 

W -- 

where 



H,W 



1 

+ 2 



H,W 



W 



I 1 

""T 2-^ / , '^pdm y^m,sPm,s Prri,s""m,! 



pd s m=-l 



(4) 



(5) 



A 



pd 



nA„ 



[rir, 



l)Ap + td-ep 



(6) 



With this choice of Ap^^, the couphng between the dP'p^ and dT-^^p^ configurations vanishes in 
the first order operator Hq, W +Tpd- Neglecting the couphng with the high-energy (P'^'^p^ 
configurations in the second order operators 



Ho^W 



W 



+ di 



Tpci,W 




Pm,s + h.C. \ , 



(7) 



^ Z_^ Z_^ A \"'m,s^m,s Pm,sPr, 



s m=-l 



pd 



(8) 



where 



A = ed-ep + AdNd - ApNp , 
we end with the effective Hamiltonian for d'^p^ configuration 



1 



H, 



eff 



Ho + E E 



f2 

pdm 



{"'m.yS^'r, 



s m=-l ^Pd 

edNd + Ud + HcF + const 



Pm,sPm,s 



(9) 



with 



hcf=j: e 



t 



^'^"^d^ d 



A 



(10) 



s m=-l '-^Pd 

In the last equality of Eq. Q we have taken advantage of the fact that every state of our 
subspace includes the non-degenerate closed p-shell of ligand for which {p^\plnsPm',s\p^) = 
6m,m' and all terms concerning the ligand become constant. We have obtained an effective 
single- ion Hamiltonian (Eq.(jni)) where the action of ligand has been reduced to the one- 
particle 'crystal field' term Hcf P0|) - Let us recall that we have assumed that 



^pdm *%^ ^ 



pd 



(11) 



Within this assumption we may neglect the terms of third and higher orders in Hejf in Eq. 
(JH). It is easy to see that the second order contributions from separate ligands simply sum 
up when we do subsequent transformations to remove hopping to first order between the 
TMI and the different ligands of the nearest surroundings. Thus, the assumption (jllj) gives 
the range of validity for the superposition model in our case. Our approach also implies that 
besides Eq. fjll|) . the characteristic energies of the terms in Hd not included in Hq in Eq. 
^ are also smaller than A^^^. At first this concerns the Coulomb energies 



15B ~ 3C < \d 



(12) 
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The main point here is that in addition to the one-particle energy difference, Ap^ contains 
also the largest Coulomb parameter A^ (see Eq. fITHj) below). 

For the chosen axially symmetric geometry, the 'crystal field' is diagonal with respect to 
the angular momentum projection m. For the general relative positions of TMI and ligands 
the HcF will have the fornt^i 

2 



HcF — 2Z zZ '^rnm' d^m.,s^rn' ,s (13) 

s m,m' =—2 

Vmm' = E {HR^)Amm'Y^'^' (^„ 0,) + h{R^)B^,n'Y:^-'^' (^„ 0,) + h,{Ri)5^^.} , (14) 



where 



m' 5v4vr^224 ^224 rj _ ( -i\m' V 47^ ^222 ^222 






t 

y 



le C^-'^ = (jiJ2"^i"^2|-^^ = mi + 7712) are the Clebsh-Gordan coefficients (see e.g. Refs. 



16| ). and the F;™ are the spherical harmonics. The coefficients 



2 



5 [ Apd 



. Mfc + 1) 
6 



t. 



2 

pdn 



A, 



pd 



(15) 



depend only on the nature of ligand and TMI and on the distance between them; the standard 
notations m = a,^ in Eq. ()15|1 correspond torn = 0, ±1 respectively. The summation in Eq. 
(fT^ goes over the ligand spherical coordinates Ri,9i,(f)i. In the spirit of the superposition 
model,— the physical and geometrical informations are separated in Eq. (fT^ . Note that in 
the general case the summation in Eq. ()13|) runs over all d-states. 

The idea to use Harrison's parametrization for the calculation of hybridization contri- 
bution to CF and to Eq. (fTK|) was first proposed in Ref. |s! from perturbative approximate 
diagonalization of the mean field one-particle part of the p — d Hamiltonian. In this scheme 
the Apd has the meaning of Hartree-Fock energies difference 

Apd,HF = £d,HF — £p,HF (16) 

ei,HF = ei + Ai {ni - 1) . (17) 



In this sense the approach of Ref. |8| is close to the 'strong' CF scheme. 

In the spirit of 'strong' CF scheme, another approach was developed in Ref. ll^ There, 
the Coulomb electron-electron interaction is taken into account after the diagonalization 
of the one-particle mean-field Hamiltonian. It is thus rewritten in terms of eigenfunctions 
of cubic CF. The applicability of Racah's parametrization of Coulomb integrals is then 
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questioned. In DMS, the 'strong' CF scheme fails. It is unable to explain the position of 
incomplete c?-shell below the Fermi level, because the mean-field neglects the configuration 
interaction. The 'weak' CF scheme is free from such difficulties and our considerations show 
the way to account for covalency in this scheme. Concluding this remark, let us mention 
that comparing Eqs. (jHl), (fT?)|) and (fTTjl we see that 



^pd 



\d,HF + Ad . (18) 



This relation explains why the TMI d-level having the Hartree-Fock mean-field energy lower 
than the ligand p-level (e.g. ed,HF = —17.77 eV for Co is lower than Sp^HF = —16.77 eV for 
oxygenS.) remains incompletely filled. As we mentioned in the Introduction, the difference 
between 'strong' and 'weak' CF schemes reflects the difference of Hartree-Fock and Heitler- 
London ways of accounting for the covalency, our consideration being close to the Heitler- 
London approach. 

B. Renormalization of Coulomb, spin-orbit and Zeeman terms 

In the phenomenological CF theory, the covalency is accounted for by introduction of 
reduction factors for Racah's parameters -B, C and for orbital angular momentum matrix 
elements. The orbital moment appears in the spin-orbit interaction Hso and in the Zeeman 
term 

Hz = fiB (^.S + l) B , (19) 

where Qs = 2.0023 is the Lande's factor and fiB the Bohr's magneton. 

The canonical transformation (Eq. (Q) changes the many-body basis of the problem. 
For the sake of consistency we should transform every additional term of the Hamiltonian 
as well as any observable. The covalency reduction factors naturally occur as a result of 
the canonical transformation of the corresponding operators. The total spin operator S 
commutes with the canonical transformation operator W (jSI) and remains unchanged. Let 
us now demonstrate the appearance of an orbital reduction factor k {L ^ kL, k < 1) for 
the spin-orbit term. We begin with the canonical transformation of annihilation operators 
and with the ligand situated at (0, 0, R) 



0"m.s ^^ (^m.s ~r 



dm,s, W 



+ 2 



dm....W 



w 
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1 ~ '^■^m ) ^m,,s — ^mPm,s ; (20) 



Pm,s 



r)'^m j VrrifS ~r '^m^m,s 



where Am = tpdm/ ^pd- The apparent similarity of Eq. (f^ with molecular-orbital expression 
should not mislead the reader. We recall that it works only in the subspace of dT'p^ and 
^n+ip5 configurations and /S.pd (jHI) depends on the number of d and p electrons n, np. 

Substituting the transformed annihilation and creation operators into the second quan- 
tization expression for an operator, we immediately obtain its transformed version. The 
spin-orbit operator acquires the form 

^SO — ~^ 2-^ "'m,s ■^m,ml ^ s,sl 0'ml,sl ~\ ^ /^ Pm,s -^^m,™/ ""s.s/ Pm/,s/ 
m,m/,s,s/ m,ml,s,s/ 



r, A^ [^ r)^mjd'm,sLm,m/^s,sr[^ r^'^mi ) dmi ,sl (21) 



^dfl y^ /^i ^ ^2 \ jt Td _ A 1 

+ ~~Z~ / , / , ^rw^ml "'m,s ^m,m,l ^ s,sl ^ml,sl i 

s,s/ m,ml =-1 



where Lm,mi is the angular momentum matrix vector for a, L = 2{d) oi L = l(p) particle; 
cr is the Pauli matrix vector. The spin-orbit coupling value ^d{p),o is the free-ion (ligand) 
one (see below). For the light ligands (e.g. oxygen) ^po *^ ^d,o, and this term is usually 
neglected. 

For the ligand situated in the point with spherical coordinates [if, 6, R) in crystallographic 
coordinates, we should first go to the local coordinate system X'Y'Z' with the Z' axis 
pointing towards the ligand. The rotation is described by three Euler angles {a = (f, (3 = 
0,7 = 0} and the operators in the two systems are related by the linear transformation 

d'm,s = H D^^^ra{a, /3, 7)c?mi,. • (22) 

The matrices Dmj,m{c(, P,^) = R^l m('^'/^'7) describe the transformation of spherical har- 
monics between two coordinate systems.— After a canonical transformation in the local 
coordinate system according to Eq. (J20|) . we perform a backward rotation and obtain 



dm,s — /, -Pmi.m(0, ~0, ~V)d'rr. 



mi,s 
mi 

A2 



dm,s + Yl ^mi,m(0, -9, -if) 



mi 



^miPmi,s r\ / J ^rm,,^^^ \'~P) "i U/^m;>., 

m2 



■ (23) 
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The summation over all ligands will give the final expression for dm.,s- Generally it is very 
complicated, but for highly symmetric surroundings, it recovers a form similar to Eq. ()20|). 
For example, for the tetrahedrally coordinated TMI we have 



4.,. = (l - l^^l - l>^i) dt,„s - ]Jlxi - Ixi Pt,„s , (24) 

"eg,s = 1 -'- "~ o TT j "'^g^^ q ^'^ Peg,s y 

where the operator d^^s annihilates the electron in the state which transforms according to 
the irreducible representation F = eg,t2g] rra = a, tt again means m = 0, ±1. Instead of 
operators Pm,s, an admixture of symmetric combinations of ligand orbitals enters the Eqs. 

(1211). 

In the ground state of tetrahedrally coordinated (i^ ion, the t2g states are filled (by 
holes) and only the non-diagonal matrix element of angular momentum {t2g L Cg) enters 
the expression of the spin-Hamiltonian parameters. It means that we may substitute 

L ^k L , (25) 

^- (-5^-^0(-^^0==(-^-^0' (-) 

in the Zeeman term ()19|) and 

Uo ^^d = Hdfi , (27) 

in the spin-orbit term. 

The matrix elements of the Coulomb interaction that were not included in i^o ©) ac- 
quire prefactors that are product of four terms (1 — |A^j. Therefore, the reduction factors 
for the Coulomb interaction may differ from those of SO and angular momentum terms. 
Nevertheless, an approximate estimate may be done by multiplying the free-ion Racah's 
parameters Bq and Cq with k^ from Eqn. ^^, i.e. 

Bo^B = k'^Bo , Co^C = k'^Co . (28) 

Such an approach is completely sufficient for the given case of ZnO:Co since the lowest ^F 
multiplet is well separated from the higher multiplets and, correspondingly, the parameters 
of the effective spin-Hamiltonian (^f-factors and zero-field splitting) are not very sensitive to 
B and C (see Section V below). 
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x/a 


y/a 


z/c 


Co 











Oi 


1 





-\ + 6 


02 


1 
2V3 


1 

2 


-\ + 6 


03 


1 

2\/3 


1 

2 


-\ + 6 


04 








1 + 6 



TABLE I: Cartesian coordinates of the C0O4 tetrahedron using the lattice parameters a = 3.2427 A, 
c = 5.1948 A and 5 = 0.0076 of the host lattice ZuO^ 

III. CRYSTAL FIELD THEORY FOR THE TETRAHEDRALLY 

COORDINATED dP ION 

The effective single-ion Hamiltonian, derived in the previous sections reads 



H, 



eff 



H, 



Coul 



+ HcF + Hso + H_ 



(29) 



The Coulomb interaction within the d-shell Hcoui will not be written down explicitly since it 
is diagonal in the many-body basis corresponding to the weak field scheme. In our approach, 
the CF parameters are connected with the local environment of the TMI (see Table IJ and 
Figure H for the ZnO:Co example). It is convenient to choose the zero of energy from the 
condition TrVm^mr = ^Y^i bo{Ri) = 0. The term Hcf in Eq. (fT^ may be split into cubic and 
trigonal parts: 

Hcf = Hcub + Hfrig , (30) 

since the ideal tetrahedron with ^ = w| and 5 = is identical to cubic symmetry. In reality, 
however, - deviates from the ideal value and 6 y^ 0. The trigonal field is described by three 
parameters. There exist several different systems of notation in the literature and we will 
use here the parameters A, v, and v' like Koidl^ and MacFarlane^^ (see also Ref. |l3|). They 
are defined as a parametrization of the crystal field matrix elements Vmm' (Eq. ()13|)) in the 
one-particle basis of the trigonal coordinate system (the z axis is the threefold axis pointing 
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FIG. 1: C0O4 tetrahedron 



towards O4): 



|a:) 



212 



2\ /T I \ 

y I ~ y-^^^i 



.^\xy) 



\z^) 



1^2 
3IX 



y^) + vfl^^) 



(31) 



\\zy) 

\v) 

1^) = -\lWy) + \ll\^y) ' 

where we have the usual real rf-basis functions constructed out of the complex basis functions 
10 = '^/1t|0) o'^ the right hand side. The three basis functions \x), \y) and \z) build up the 
t2g representation of the cubic group and |f ), \w) span up the Cg subspace. The one-particle 
crystal field matrix elements are given in this basis by: 



zz 


— 


5^- 


3 




r 

XX 


= 


T/ — 

^yy - 


iA + 


V 

3 


r 
vv 


= 


v^^ = 


--1^ 




r 

XV 


=^ 


v' . 







(32) 



The relationship between the A, v, and v' parameters with the Stevens equivalent operators 
and other parametrizations used in the literature can be found in the Appendix. 

If we substitute the ligand coordinates from Table E into Vmm/ in Eq- ffTH) and then 
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transform them into the trigonal crystal field parameters we obtain 



A = -^ Mi?i 



72(^2^2 _ ^4) _ 3^4 ^ 20y62i( 



24:Rt 



V 



vt 



9Rf 
3 

~14 

5b4{Ri 
9Ri 

V2 
14 



20 



3{zfa' - zt) 



5v/6, 



Zia^ 



hiRi 
20 



63 



64(i?4 



b2{Ri)^^^ + 2b2{R,) 



V2 
7 



3(zi'a2 - 4) 



(33) 



Vs. 



2;ia 



x/2 



12 



64(^4 



+ ■ 



&2(i?l)^^^^ + 262(i?4) 



where Zi is the z coordinate of the ligand Oj; -Ri = Ja'^/S + zf, R4 = Z4 are the correspond- 
ing distances. 

The importance of the non- diagonal matrix element between Cg and t2g states, v', was 



Q 



first pointed out in Refs. |1^ and |20| and will be outlined in the following, but it was not 
thoroughly treated in the standard text books.— Hopping integrals tpdm have been calculated 
from Harrison's table. ^ Then all the CF parameters depend only on one value: Apd- Its 
determination is complicated by the fact that the Coulomb repulsion Ad is partially screened 
in semiconductors and it differs much more from free ion value than B and C. That is the 
reason why we have used Apd as an adjustable parameter. After determination of its value 
from experimental knowledge of cubic parameter A, the other parameters v and v', are 
determined by geometry. 

For the spin-orbit and Zeeman terms we have 



Hso 
Hz 



r) / y m,s 



Hb [9sS + kVj B 



J-^m.mJ tr Q 0/ Cin 



^m.ml ^ s.s/ ^mLsf 



(34) 
(35) 



here ^d is the renormalized spin-orbit coupling, and k is approximately given by Eq. (j26|) : S 
and L are respectively the total spin and orbital angular momentum operators. Within the 
^F term, Hso ^lay be rewritten as 



Hso,^F = ASL , with A = — ^d/3 . 



(36) 
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IV. SPIN HAMILTONIAN FOR THE GROUND STATE MANIFOLD 

As it was already noted in the Introduction, this second step of our calculation can be 
performed exactly by numerical diagonalization of the effective single-ion Hamiltonian ()29|l 
once we have determined all its parameters (see Appendix). In that sense it is not restricted 
to the 'weak' field case. In addition to the numerical diagonalization, we derive below an 
analytic perturbative formula which is valid in the 'weak' field case of ZnO:Co. 

A. Many-body basis 



For a while let us neglect the spin-orbit coupling. Then the total spin S and total 
angular momentum L are conserved, because Coulomb interaction is rotationally invariant 
and does not depend on spin. The ground state for seven rf-electrons (three holes) is \^F) 
with L = 3, S = 3/2 and it is (2L -|- 1)(2S' -|- 1) = 28-fold degenerate. The eigenfunctions 
are |L,M, m^), M being the momentum projection on z axis, rris the total spin projection. 
For M = 3 and m^ 



I we have 



3. 3, 1 



7t ^t ^t 



Cinl-Cil -t-Ciot IVCXCj 



Acting on this state by I/_ operator we obtain successively 

3,1, 



3,0,- 




'^-2t'^lT 2t "I" \ p."-lf"0T'^2t I 



\vac) 



I ZiLL o-t-Ciri-t- ttoi- ~n a 1-1- ^OT IT/ t/CiC) . ot-C 



Under the action of cubic crystal field this level splits into 1 singlet and 2 triplets. Then the 
basis functions may be labeled as |3,x, ^, '^s), where x = ^2,^2,7"! denotes the represen- 
tation of the cubic group, m is the projection of a fictive angular momentum within each 
manifold. We have 

13,^2, 0,772^) = — — |3,0,m,) + — (|3,3,m5) - |3, -3,m,)) 

1 [5 

|3,T2,±l,m,) = — =|3, =F2,m^) ± W-|3,±l,m,) 
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\3,T2,0,ms) = -y={\3,3,ms) + \3,-3,ms)) (37) 

/5 1 

|3,Ti,±l,ms) = =Fy-|3,=F2,m^) + ^ |3,±l,ms) 

2 1/5" 

|3,Ti,0,m,) = --|3,0,m,) --W-(|3,3,m,) - |3, -3,ms)) . 

The trigonal field splits the triplets into doublets and singlets and also couples the states of 
different manifolds with equal m. This is schematically shown in Fig. |21 The lowest excited 
level is ^P with L = 1, ^ = 3/2 



1,1.5 

' '2 



J -dl2'\d{^dl^ - Y -c?Li^4t4t I^«c) 



1, 0, -\ = W - (^dl^idl^dl^ - 2dl^^dl-^d\^'^ \vac) , etc. 

In the free ion it is separated by energy 155 from the ground state. The cubic field couples 
the states \l,M,ms) and |3,Ti,m, m^) with m = M. The trigonal field has also matrix 
elements between |3, A2, rh, rris), |3, T2, m, m^) and |1, m, rrts) states that are proportional to 
V and f', thus they are much smaller. 

B. Perturbation theory 

From the parameter values discussed below, it will become clear that we are in the regime 
where 155 ^ A ^ v,v',\. In the strong cubic crystal field case (i.e. when A ^ 155) a 
perturbative formula for zero field splitting 2D and the gyromagnetic factors g\\ and g^ was 
developed by MacFarlane .^^'^^ But it is not adequate in the present situation. Now, we can 
define four small values 7^5^, K ~ Z ~ K" ^"^ fact, we have 155 ^ A, and the value A/155 
being almost of the same order of magnitude as f /A, then the ratio 

V V A 

155 ~ AT55 

can be considered as an order of magnitude smaller than A/155. This justifies the applica- 
tion of the weak crystal field approach. 

We will proceed in three steps that may be regarded as three subsequent canonical trans- 
formations of our Hamiltonian similar to Eq. (JH). First, we eliminate the coupling with ^P 
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T, 



'E 



! E, 



'1/2 



'1/2 
'3/2 



E 



3/2 



'1/2 



F 



E 



A, 



'1/2 
'1/2 



'1/2 



HE 



3/2 



'1/2 



-3/2 



Free Ion 



Td 



C3y 



'1/2 



SO 



FIG. 2: Qualitative splitting due to the cubic field, followed by the trigonal crystal field and by 
the SO effect. 



states retaining only the order A/15-B (the explicit use of weak CF scheme). Then only the 
states |3,Ti,m, -m^) acquire an admixture 



\Ti,m,ms) = \3,Ti,m,ms) + 



(1, m, rusl HcF |3, Ti, m, m^) 



where 



(-^TiO ~ Epo) 

l,m, msl iJcF |3,Ti, 771,7713) 2 A 



l,m,ms) , (38) 



= Kp 



(Et.o - Epo) 5 15B 

We thus obtain an effective Hamiltonian acting in the ^F subspace. In the next step we 
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consider the perturbation due to the trigonal field with the small parameters -^ and ^. In 
the following we will use the first order ground state wave function of the approximately 
diagonal crystal field Hamiltonian,— i.e. 

\iPq) = \A2, 0, nis) ^ |3, A2, 0, nis) - kq l^^i, 0, m^) , (39) 

as well as the excited states 

|T2,0,m^) ^ |3,r2,0,m,) 

|T2,±l,m,) ^ \3,T2,±l,m,)T'^2\Ti,±l,m,) (40) 

\Ti,0,ms)t ^ |Ti,0,ms) + Ko 13,^2, 0,ms) 

|Ti,±l,m^)j ^ \Ti,±l,ms) ± K2\3,T2,±l,ms) , 

where 



^Vio ,, ^ , t^Vio,, ^ , 

Ko = —^ ^— (1 + 2«:p) ~ — — (1 + 2^p) 

5 {Et^o - Ea^) 9 a 

v + 2v'^ Up ( 3^2 A 

K2 = — 1= 1 1= V V 

2Vh{ET,i-ET,i) V^{Et,i-Et,i)\ 2 ) 

The corresponding energies are 

Et,i ^ -^ - ^ (41) 

-EtiO 
-E'tiI 

Then we consider the spin-orbit interaction as a perturbation with respect to the crystal 
field Hamiltonian. We thus obtain the usual formulae for the (yf-factor and the anisotropy D 

9f.u -9s = -2XkA^^ , (42) 

D = -X^ (A,, - A,,) , (43) 



5 




A V 




1 — 




5 3 




A V 




5 6 




3A 3v Av'V2 


4 A2 


\ \ 


— 


5 5 5 


5 755 


3 A 3v v'V2 


4 A2 


+ 




5 10 210 


5 75fi 



where 
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with Eq = Ea^ and A is defined in Eq. fl36|) . These are the parameters appearing in the 
effective spin (5* = |) Hamiltonian: 



Hspin = fJ'B9\\BzSz + I^bQi. (B^Sx + BySy) + D 



S^^--S{S + l) 



(45) 



Let us note that all energy denominators appearing in A 



/lU 



are of the order of cubic 



splitting A. Thus, the perturbation theory requires the SO coupling to fulfill A ^ A. But 
A may be of the same order of magnitude as v and v'. The operator L^ couples the ground 
state only with |T2, 0, m^) 



L, 1^2, 0, m,) = (2 + KoVs) IT2, 0, m,) 



then 



For L, 



A, 



4 (1 + K0V5 



{E_ 



T20 



E 



A2) 



{U 



-Ko 



^ LS)/2 we have 

L, \A2, 0, m,) = - {-2V2 (IT2, 1, rris) - IT2, -1, m,)) 

--^ (|Ti, 1, m,) + |Ti, -1, rus)) -J- (IT2, 1, m,) - jTa, -1, m,)) 



and 



A ~ 



4 1 



2 



The anisotropy constant appears only in the third order (second order of spin-orbit and first 
order of trigonal field) 



D 



AX' 



V 



A V2A 



3kov5^ 



A2 

A 



'a 



10^2 t;' 
3 A 



1 + 2kp) 



The final results are: 



D = ^ 



^A 3 A 



75 B 



9\\-9s ^«;[i 3A+9A^V^ + 75i? 



(46) 



91. = 9s 



8A 



A; 



-^ ^ 6A ISA*^ V 75B 



An alternative perturbative formula for zero-field splitting 2D^ which is valid in the present 
situation, was derived by Mao-Lu and Min-Guang.— However, our result is much more com- 
pact than theirs (the Eqs. (5)- (9) of Ref. llJtake one page and a half) and, correspondingly, 
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0,03 



0,02 



0,01 



0,00 



perturbation theory 

MacFarlane 

exact diagonalization 




A/B 



FIG. 3: Comparison between our perturbation theory (|46|) for the weak field case (dashed hne), 
the result derived by MacFarlanei^ in the strong field case (dotted), and the exact diagonalization 
calculation for Sg = g± — gy (full line). 



more practicable. We have checked that the difference between the compact formulas ()4fj|l 
and the result given in Ref. |lj is very small and can be neglected in numerical applications. 
To check the applicability of our perturbation theory PUJ) we compared it with the exact 
diagonalization (see Appendix). We fixed the parameters (in inverse centimeters) v = 54, 
v' = -213, C = 3148, B = 804, k = 0.85, and k^d,o = 481 to those values which we derived 
for Co^+ in ZnO (see next Section) and varied A/B to display the different regimes. As one 
can see from Fig. El the difference of ^f-factors Sg = g± — g\\ from our formula compares well 
with the exact diagonalization for A/B up to 15. Therefore, we can certainly apply it to the 
case of Co^"^ in ZnO where A/B is roughly 5 and where there are considerable deviations 
for MacFarlanes formula.— On the other hand, the exact diagonalization converges towards 
MacFarlanes formula for large values of A/B. 
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V. APPLICATION TO ZNO:CO 



We apply our theory to the calculation of spin-Hamiltonian parameters for Co impurity 
in zinc oxide. As input, we need the following parameters of H (Eq. (^): (i) the structure of 
the Co environment given in Table HJ (ii) the charge transfer energy ApwfEq. ()18|l ). (iii) the 



^iv) the Racah's 

Ii 



free-ion spin-orbit coupling ^^q ~ 567 cm^^ (see e.g. Table 7.6 of Ref. |3|)j 
parameters B^ ^ 1115 cm~^ and Cq ~ 4366 cm~^ (from Table 7.5 of Ref. '3). 

As was mentioned above, we calculate the hopping integrals that enter H (Eqs. (^ and 
(jHI) from Harrison's expressionsi^ 

tpdm{R) = Vpdm P7/2 ' ^P'^'^ = -2.95, T^p^^ = 1.36, (47) 

where the value r^ = 0.76 A for Co ion. The distance R is measured in A and tpdm in 
eV (1 eV = 8065.5 cm^^). This gives, e.g. for tpda^Ri) ~ 1.34 eV. The coefficients 
bk (Eq- (HSI)) are inversely proportional to the charge transfer energy Ap^. We choose 
the Apd ^ 3.6 eV ~ 28800 cm^^ so that the cubic splitting A (Eq. (J33|) ) is equal to 
the experimentally determined valued A = 4070 cm~^. Then, the trigonal CF parameters 
t; ^ 54 cm~^ and v' ~ —213 cm~^ are unambiguously determined by the Co environment 
via Eqs. ()33p and are not additionally adjusted. The SO and angular momentum reduction 
factor k ^ 0.85 was calculated iiom. Eq. (jlHI)- This gives the spin orbit coupling ^^ = kC,d,o ~ 
481 cm~^, A = — ^d/3 ~ —160 cm^^ very close to the values met in the literature (and 
adjusted empirically): ^d = 450 cm~^ in Ref. tlJ, and ^d = 430 cm~^ in Ref. |2. From our 
above consideration, we have seen that the reduction factors for Coulomb interaction may 
differ from those of SO and angular momentum terms. Nevertheless, an order of magnitude 
estimate may be done as i? ~ k'^B^ ^ 804 cm~^ and C ~ k'^Co ~ 3148 cm~^. We adopt 
these values for our calculations. The agreement with the experimentally adjusted values 
B = 750 cm~^, C = 3500 cm~^ of MacFarlanoi^ is very good, keeping in mind the roughness 
of the estimate. We should also note that the influence of the B and C parameters on the 
calculated g-iactors and the zero field splitting 2D is rather small. So, using instead of our 
estimations for B and C the values of MacFarlane (and keeping all the other parameters 
constant) leads to 2D = 4.882, g^\ = 2.230, and g± = 2.256, quite close to the results listed 
in Table II. 

The parameters of the ground state spin-Hamiltonian (Eq. ()45|) ) obtained by exact di- 
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P. Koidl^ 


R.M. MacFarlanei^ 


Present work 


experiment 


V 


-120 


-400 


54 




v' 


-320 


-350 


-213 




A 


4000 


4000 


4070 


4070^ 


B 


760 


750 


804 




C 


3500 


3500 


3148 




k 


1 


0.8 


0.85 




Hd,o 


430 


450 


481 




2D 


5.44 


5.41 


4.91 


5.52" 


9\\ 


2.24 


2.20 


2.23 


2.236ii 


9± 


2.28 


2.23 


2.26 


2.277iii 


2D' 


21.4 


90.3 


5.8 


382 


9\\ 


2.85 


3.36 


4.22 


3.5223 


^■2EE - E^AE 


15171 


15050 


14381 


151232a 



TABLE II: Measured ESR^i and magneto-optiG23 data, compared to those calculated from CF 
theory with empirical parameter32frW or estimated from our approach. 2D {2D') and g {g') are the 
values for the ground (excited) state. The energy unit is inverse centimeter. In the empirical CF 
theorySii^ all the 7 parameters listed in the upper part of the table are used as fitting parameters, 
whereas they are microscopically calculated in our approach fixing only the cubic CF splitting A 
to the experimental value. 



agonalization of single-ion Hamiltonian Heff (Eq. ((221)) ai'e shown in the Table ITTl Note 
that our sign convention for the crystal field parameters corresponds to the electron repre- 



I2I . The agreement with 



sentation in contrast to the hole representation used in Refs. |2| and 

purely empirical approaches, where all parameters are fitted to experiment, is very good. 

We show also the parameters for the excited state. 

Table IIIII compares the spin-Hamiltonian parameters obtained by analytic perturbative 
approach (Eqs. ()46p) with the results of the exact diagonalization and experiment. We see 
that the analytic results lie within 15% of accuracy. For the phenomenological parameter set 
of Koidl^ the accuracy is about 20%.— The reason is that in our set the absolute values of 
trigonal parameters are smaller than for the set of Ref. Ij and the cubic splitting is the same. 
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experiment 


diagonalization 


perturbation theory 


2D 


5.52 


4.91 


4.31 


9\\ 


2.24 


2.23 


2.25 


g± 


2.28 


2.26 


2.28 



TABLE III: ESR data as calculated from the numerical diagonahzation and the perturbation 
theory using the calculated crystal field parameters. 



thus the perturbation theory for our set converges better. The main reason of the deviation 
from exact diagonahzation is our neglect of the interaction with ^G term that violates the 
Hund's rule, but nevertheless lies rather low in energy, just above the ^P term (see Table 
rvj) . This deviation is remarkable for 2D, but much less for the (^-factors. 



VI. CONCLUSION 



We have shown that Harrison's parametrization of electronic structure of solids^ may be 
successfully applied to the calculation of spin-Hamiltonian parameters (Eq. (03))) for TMI 
impurities in semiconductors (Eq. ()46|)). It is especially useful for the description of low 
symmetry paramagnetic centers as it provides the connection between CF parameters and 
the geometry of TMI surroundings (Eq. (jHSI))- Thus, the number of empirically adjustable 
parameters is substantially reduced. 

We have demonstrated that the physical reason for the possibility to apply the CF concept 
to TMI in semiconductors is the strong Coulomb repulsion within the d-shell. It provides 
the large value of charge transfer energy Ap^ (Eqs. © and (jlSp ) even in the case when the 
mean-field energy of the d-level falls into the valence band. We have given the explicit form 
of the canonical transformation (Eqs. (^ and ©) of the many-body Hamiltonian (Eq.(^) 
and basis functions, which exploits this strongly correlated feature of the TMI subsystem 
(Eq. (fTTj) ). and provides the effective single- ion Hamiltonian (Eq. Q). The latter connects 
the CF Hamiltonian with the geometry of local surroundings of the impurity and with 
the parameters of the electronic structure (Eqs. p3p to ()15|)). The transformation (Eqs. 
(0} and ((HI)) accounts for the covalency in the 'weak' CF case within the Heitler-London 
configuration interaction approach. When applied to the spin-orbit, Zeeman and Coulomb 
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terms, it renorinalizes their parameters by covalency. 

We have apphed this theory to the Co impurity in ZnO. We have adjusted only one 
parameter of our starting p-d Hamiltonian (Eq. (0)), which acts in an energy scale of several 
eV. In the result, we have fairly well reproduced a number of measurable quantities available 
from ESR and optical experiments. Note that these values reflect the tiny features of 
electronic structure (magnetic anisotropy, Zeeman splitting), which have the scale of several 
cm~^. The results indicate that the proposed theory catches the essential physics of TMI in 
semiconductors. 
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Appendix: Exact diagonalization for 3d^ (3(i^) 

In this section we detail the numerical calculations of the exact diagonalization for 3 
particles on a 3d-level. The case of interest namely Sd"^ can be obtained from the former 
using a particle-hole transformation. 

The number of states for 3 particles on a c? level is 10 x 9 x 8/3! = 120. These states can be 
labelled by usual quantum numbers for total spin and angular momentum, we thus obtain 
different multiplets: "^F for which S" = |, L = 3, this multiplet contains (2S'+1)(2L + 1) = 28 
different states. We also have ^P {S = ^,L = 1),^H {S = ^,L = 5),^G {S = ^,1 = A), ^F 
(5 = i, L = 3), 2/} and ^D' (S = |, L = 2), and finally ^P (5 = 1, L = 1). This basis will 
be noted {\a) = \S,L, M, rris, e)}, the last quantum number e is needed to distinguish states 
belonging to the 2 multiplets ^D and ^D'. 

The {\<y)} basis is the natural one for the Coulomb interaction as well as the Zeeman 
Hamiltonian, however an other basis emerges when writing the one-body part of the Hamil- 
tonian, namely the crystal-field and spin-orbit terms. Let us denote by {\i),i = 1..10}, the 
basis for one electron on a. d level, i is an index for the (mi, a) state, where mi and a are 
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Term 


Coulomb energy 


^F 


3A - 15B 


4p 


3A 


^G 


3A-11B + 3C 


^H 


3A-6B + 3C 


2p 


3A-6B + 3C 


2p 


3A + 9B + 3C 


^D, 2D' 


3A + 5B + 5C± v/(193S2 + 8BC + AC^) 



TABLE IV: Coulomb energy depending on the multiplet 



the momentum and spin quantum numbers. We can construct a new basis of 120 states 



leveL 



JJJ 



ijk) = CiCjcll vac), i < j < k : 1..10, n : 1..120} where | vac) is the empty d 



The complete Hamiltonian has been written in Eq. (j^Uj) . The Coulomb part is diagonal in 
the {\a)} basis, and the Tab. IIVI gives the different energy values. Contrary to MacFarlane's 
work^i where the Zeeman Hamiltonian is treated in a perturbative manner, it is here diag- 
onalized on the same foot as the other terms. In the {\a)} basis, the Zeeman Hamiltonian 
is block-diagonal, and the non zero matrix elements just connect states by !'+(-) or 5'_|_(_). 
The crystal field Hamiltonian is the sum of the cubic and trigonal parts. The one particle 
matrix elements of Hcf (Eq. CHI)) may also be expressed in terms of Stevens equivalent 
operators^ 

corresponding to a d} configuration. The Stevens operators are given by 



--fi0(6°-20V2O^) Ht„g = B'^dl + B'^d\ 



(48) 



01 
01 
09, 



35Ll - 30L(L + l)Ll + 25Ll - 6L{L + 1) + 3L\L + 1)^ 

1 

4 

3Ll - L{L + 1) 



■{Lm+L^_) + {Ll^Lt)L,] 



(49) 



where the operators Lz,L^, or L_ are one particle operators. 

In analytic calculations we have used the {|«)} basis and for the rf^ configuration the 
crystal field Hamiltonians (Eq. ()48|)) have to be used with the parameters 



-D. — :—. 



b: 



b: 



B' 



B' 



(50) 
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The parameters A, v, and v', previously defined in Eq. ()32p. are connected with the Stevens 
parameters by: 

-°4 - 120 360 '^^^ 2 '^ ) 

B', = -U^+'-^v') (51) 



140^ 

o/ _ t;-2V2t;' 
-°2 — 21 



Using Eq. (jHHj) we obtain 



Bl = -^b,{R,f-4 
^ 2880 ^ ^ i?t 

R, 1 i^^^^ 7V6z,a' + 72{zt-zy)+3a' , &4(i^4) ] .„. 

^^ = "140 1^^^^^^ 216^^ + ^1 ^^^^ 



2lo{^^(-^)^"^ 



B2 = -;^^2(i^l)^^5^ + 262(i^4) , 



where Zi is the 2; coordinate of the hgand 0^; -Ri = Ja'^/S + zf, R4 = Z4 are the correspond- 
ing distances. For completeness we give here also the relation with another parameter set, 
that is often met in the literature 

Dq = 125°, Dt = 12BI Da = -SB'^. (53) 

This set is used e.g. in Ref. Il^ 

The crystal field Hamiltonian is easily written in the one-particle basis {\i),i = 1..10}, 
where one evaluates the matrix elements Hf-^ = {i\HcF\j)- Then the matrix elements of 
the CF Hamiltonian can be written in the 3-particle basis {\n) = \ijk)} as follows 

{k'j'i'\HcF\ijk) = H-/- {6j'j6k'k — Sj'kh'j) — H-/j {6j'i6k'k — Sj'kh'i) + H^^, {6j'i5k'j — Sj'jh'i) 

~B^j'i i^i'j^k'k — ^i'k^k'j) + Hjij {6in5k'k — ^i'k^k'i) — Hj,f. {6i'i6k'j — ^i'j^k'i) 

+-^fc'i i^i'j^j'k — ^i'k^j'j) — Hi^,j {6i'i6j'k — Si'kSj'i) + Hj^,}. {6i'i6j'j — 6i'j6jfi) 

where 6 is the Kronecker symbol. The spin-orbit term (Eq. (jS3)) is also easily written in the 
one-particle basis {\i)}, then in the 3-particle one {\ijk)}, using the preceding expansion. 

At this stage we have some part of the total Hamiltonian written in the {\a)} basis, the 
other one in the {\ijk)} one. To perform the numerical diagonalization, the last quantity 
needed is the transformation matrix to connect these two basis. The transformation basis is 
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a kind of Clebsh-Gordan coefficient matrix for three particles constrained by Pauli principle. 
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